Enhancing quantum entanglement by photon addition and subtraction 
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The non-Gaussian operations effected by adding or subtracting a photon on the entangled op- 
tical beams emerging from a parametric down- conversion process have been suggested to enhance 
entanglement. Heralded photon addition or subtraction is, as a matter of fact, at the heart of 
continuous-variable entanglement distillation. The use of such processes has recently been experi- 
mentally demonstrated in the context of the generation of optical coherent-state superpositions or 
the verification of the canonical commutation relations. Here, we carry out a systematic study of 
the effect of local photon additions or subtractions on a two-mode squeezed vacuum state, showing 
that the entanglement generally increases with the number of such operations. This is analytically 
proven when additions or subtractions are restricted to one mode only, while we observe that the 
highest entanglement is achieved when these operations are equally shared between the two modes. 
We also note that adding photons typically provides a stronger entanglement enhancement than 
subtracting photons, while photon subtraction performs better in terms of energy efficiency. Fur- 
thermore, we analyze the interplay between entanglement and non-Gaussianity, showing that it is 
more subtle than previously expected. 

PACS numbers: 03.67.-a, 03.67.Bg, 42.50.-p, 42.65.Lm 



I. INTRODUCTION 



Quantum information processing with Gaussian con- 
tinuous variables is a well-established subfield of quan- 
tum information sciences today, see e.g. pJ-[|J]. Quantum 
key distribution, for example, can be carried out dealing 
with Gaussian states and measurements only . Nev- 
ertheless, non-Gaussian quantum states and operations 
arc indispensable to perform certain other continuous- 
variable quantum information tasks, such as quantum 
entanglement distillation , quantum error correction 
or universal quantum computation (Tlj . In addi- 
tion, any Bell test of quantum nonlocality that relies 
on Gaussian measurements only necessarily requires the 
preparation of a non-Gaussian entangled state p^4l6| . 
while a quantum bit commitment protocol that is se- 
cure against Gaussian cheatin g m ust necessarily involve 
a non-Gaussian resource state [T7| . 

Deterministically producing a non-Gaussian quantum 
optical state by using the Kerr effect is unfortunately un- 
feasible today because it requires quite high (called giant) 
optical nonlincarities, which are not accessible in the lab- 
oratory. Probabilistic degaussification schemes, however, 
have been shown feasible based on photon addition and 
subtraction. Not being unitary, photon addition or sub- 
traction can only be achieved probabilistically, that is, 
conditioned on a particular measurement outcome. One 
thus refers to it as heralded photon addition or subtrac- 



tion. The effect of photon subtraction can be obtained 
by sending a small fraction of the optical beam on an 
avalanche photodiode and conditioning the use of the re- 
maining fraction of beam upon photon-counting events 
fTsL . Photon addition can be achieved as the result 
of a single-photon excitation of the light field produced 
by parametric down-conversion in a nonlinear medium, 
conditioning the use of the signal output mode on the 
detection of a photon in the idler mode [2(j . 

In principle, an arbitrary single-mode state can be pre- 
pared by applying a sequence of photon additions [22[ 
or subtractions [23[ that are interleaved with displace- 
ment operations; similarly, an arbitrary operation de- 
pendending only on the photon number operator can be 
generated by using appropriate superpositions of addi- 
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tion and subtraction [2 If . On the experimental side, 
the use of photon subtraction from a squeezed vac- 
uum state has been demonstrated in [2J-[2fJ| in order 
to generate low-amplitude coherent-state superpositions 
(Schrddinger kitten states) of traveling light (schemes 
that can be further developed to generate larger-size 
squeezed Schrddinger cat states [30|), while the use of 
photon addition combined with displacements has al- 
lowed to generate arbitrary superpositions of the first 
three Fock states [3l|. Moreover, the ability to su per pose 
different sequences of additions and subtractions j32[ has 
enabled checking of the canonical commutation relations 

mm- 

In this paper, we focus on the enhancement of quantum 
entanglement that results from adding or subtracting an 
arbitrary number of photons on the two beams emerg- 
ing from a nondegenerate parametric down-conversion 
process, in an attempt to understand the generally 
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FIG. 1: (Color online) Degaussification schemes analyzed in this work. Alice and Bob share a two-mode squeezed vacuum 
state, which results from applying a two-mode squeezing operation S(r) on two vacuum modes (A and B). They perform k 

and I photon additions (a) or subtractions (b) on their corresponding modes, generating the output states IV'idiP) or IV'sub'); 
respectively. 



admitted — but not systematically analyzed — notion that 
degaussifying the down-converted beams makes them 
more entangled fiH, I35rl43j . In Section II, we provide the 
basic equations describing the state obtained by adding 
k and I photons on the first and second beam of the two- 
mode squeezed vacuum state, respectively, or when sim- 
ilarly subtracting photons. Section III is focused on the 
case I = 0, where we can analytically prove that quan- 
tum entanglement is a monotonically increasing function 
of k. Section IV treats the general case (k,l > 0), and 
presents an exhaustive analysis of the behavior of entan- 
glement enhancement as a function of k and I. Note that 
a subclass of these states (k = I) has already been ana- 
lyzed in (U H3|. In Sec. V, we investigate a measure of 
the non-Gaussianity of these photon-added and photon- 
subtractcd states, and show that the link between this 
measure and entanglement is rather subtle. Section VI 
is devoted to our conclusions. 

We would finally like to remark that, even though we 
are focussing on addition and subtraction of photons in 
an optical mode, our analysis applies also to addition 
and subtraction of excitations in other platforms such 
as phonons in mechanical oscillators [44[, or polaritons 
in atomic ensembles [45j . In particular, atomic ensem- 
bles might offer several advantages over photonic sys- 
tems pa. l47j: first, they act as a quantum memory, and 
hence, the state can be re-used if the addition or subtrac- 
tion protocol does not succeed; secondly, generating the 
two-mode squeezing interaction between an optical mode 
and the ensemble is not more challenging from the exper- 
imental point of view than generating the beam splitter 
interaction, which means that addition and subtraction 
are on an equal footing in the laboratory. 



II. ADDING OR SUBTRACTING PHOTONS 
ON A TWO-MODE SQUEEZED VACUUM STATE 

Our starting point is the two-mode squeezed vacuum 
(TMSV) state. Let us call A and B (for Alice and Bob) 
the involved modes, whose boson operators are respec- 



tively denoted by {a, a'} and {b, b'}. These boson oper- 
ators satisfy the usual canonical commutation relations 
[a, a*] = [6, 6+] = 1. The TMSV state is obtained by 
applying the joint operation S(r) = exp[r(a& — a^w)/2] 
with r G [0,+oo[, known as the two-mode squeezer, to 
the vacuum state of modes A and B, that is, 



|TMSV) = S(r)|0,0) = yjl - A 2 ^ A r >,n), (1) 

where A = tanhr G [0,1[, and |m,n) = \rn)A <& \ti)b- 
Here, {|n)} nS N denotes the number states defined by 
a\n)A = \fn\n — 1)a and at|n)ji = y/n + l\n + 1)a for 
mode A, and analogous expressions for mode B. 

Our aim now is to study how photon addition and sub- 
traction affects the TMSV state when applied locally by 
Alice and Bob. To this end, we consider the schemes de- 
picted in Fig. [T] In the first scheme (Fig. [T^,), Alice and 
Bob add, respectively, k and I photons to their mode. It 
is straightforward to show that the final, properly nor- 
malized state can be written as 
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where 2-Fi(a, b; c; z) is the Gauss hypergeometric function 
defined as a series expansion 

_ , , . ab a(a + i)b(b + l) , , A . 

2 F l{ a, b ;c;z) = l + -z + ^ + . . . (4 ) 

In the second scheme (Fig. [TJd), photon addition is re- 
placed by photon subtraction, and the output state can 
be written as 



i^b > = Ev«^'V-*> n -*>, (5) 
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2 F l (k + l,k + l;l + k-l;\ 2 ) X (*) 



(6) 



where we have assumed that k > I (exactly the same 
expression but interchanging k and / holds for k < I). 

Note that here we treat photon addition and sub- 
traction as ideal (a\a) operations. In realistic schemes 
based on the beam splitter (for subtraction) or the two- 
mode squeezer (for addition) interaction with an auxil- 
iary vacuum mode, this is an approximation which be- 
comes exact only in the unphysical limit of vanishing 
interaction (e.g, for perfect transmissivity of the beam 
splitter). Nevertheless, as long as the interaction is 
kept very weak — which then makes successful subtrac- 
tion or addition events rare, but still frequent enough 
for applications — the idealized description is a good ap- 
proximation |32( . In any case, we refer the reader to 
[2ll . l23l . [30l |42| for a rigorous treatment of photon ad- 
dition and subtraction in experimentally realistic condi- 
tions. 

In the following, we analyze the entanglement of these 
states as a function of the number of photon additions 
or subtractions. Being pure bipartite states, their entan- 
glement is uniquely measured by the entanglement en- 
tropy [iq . defined for an arbitrary state as = 
5[trs{|?/')(i/)|}], where S[p] = — tr{plog/5} denotes the 
usual von Neumann entropy of the density operator p. 
In our case, evaluating this quantity is straightforward 
because Eqs. ^ and §5§ are in Schmidt form, so that the 

entanglement entropy of the states IV'idcP) an d l^iub) * s 
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and 



k,i) 



(7) 



= - E ff^bgg^, (8) 

respectively. Unfortunately, we have not been able to 
carry out these sums analytically except in the trivial 
case k = I = 0, where we get the well-known expression 
for the entanglement entropy of the TMSV state [|| , 



TMSV 



(A) 



1-A 2 



los 



A 2 



log 



1 - A 2 



(9) 



which is a monotonically increasing function of A. De- 
spite this absence of a closed expression for Eqs. (|7|) and 
([5]), we are able to analytically derive their dependence 
on the parameter k for I = in the next section. 



III. ENTANGLEMENT PROPERTIES FOR 
ONE-MODE OPERATIONS 

Let us restrict to the case in which only one of the 
modes undergoes photon addition or subtraction opera- 
tions, while the other is unchanged. In that case, we will 




FIG. 2: (Color online) Three equivalent degaussification 
schemes. Starting with the vacuum state, each of these three 
single-mode operations results in the same output state 
up to a normalization factor (i.e., the probabilities of success 
are different): Alice adding k photons after the two- mode 
squeezer S(r) (solid-border pink box), Bob subtracting k pho- 
tons after S(r) (dashed-border green box), or Alice adding k 
photons before S(r) (yellow box with dashed-dotted border). 



be able to prove analytically that the entanglement en- 
tropy, either £^dd^ or ^iub^' can on ly increase with the 
number of operations k. The first thing to note is that 
the three schemes shown in Fig. [5] lead to the exact same 
state 



\ip {k) ) = J2\/p { n k) \n + k,n), 



with 



P 



« = (1 - A 2 ) fe+1 A 2 



n + k 



(10) 



(11) 



Note that this state follows from by putting I = and 
using 2-Fi(fe + 1, 1; 1; A 2 ) = (1 — A 2 ) - * -1 . Using Eq. Q, 
it is easy to prove that 



6|TMSV) = AatlTMSV), 



(12) 



implying that Alice adding k photons on the first mode 
(solid-border pink box in Fig. [2]) has the same effect as 
Bobs subtracting them from the second mode (dashed- 
border green box in Fig. [2]) , up to a normalization factor 
related to the success probability of the corresponding 
operation. Secondly, it is also easy to check that 



S(r) (at) fc |0,0) 



1 



(at) S(r)|0,0>, (13) 
cosh r 

where we have used Eq. (fT2")) as well as the relation 
S'(—r)a'S(—r) = a^coshr — 6sinhr. Hence, adding k 
photons before or after the two-mode squeezer is equiv- 
alent, except for a normalization factor again related to 
the probability of success of the non-unitary operation. 

Since Eq. (fT0|) is in the Schmidt form, the entanglement 
entropy of state |?/>( fe )) is easily evaluated as 



E (k) 



(fe) 



(14) 



71=0 



4 



V A = 0.4 



(a) 



\\ < 



3.5 



X/ 

/// 



k 



10 




A; 



10 



4.0 
3.5 
3.0 
2.5 
2.0 
1.5 
1.0 



4.0 
3,8 
3.6 
3.4 
3.2 
3.0 
2.8 
2.6 



* 

* 

* 

* 

i 

f 


(c)" 


2 4 


6 8 10 


5 / / 
o / 

/ 


(d) 

\ X 

*s 

\ \ 
\ \ 


5 , / 

i 

/ / 
* y 

3 1 ( 


\ \ 
\ \ 
\ \ 
\ \ 
"v \ 

\ 


2 4 


6 8 10 


:JS?(M) 

2 4 




•;.*, • 

k 

6 8 10 



FIG. 3: (Color online) Entanglement entropy of the photon-added and -subtracted states, I'^dd' 1 ) an d IVsub )> as a mnc tion 
of the number of operations (k, I) for A = 0.4. Note that we have normalized it to the entanglement of the corresponding 



TMSV state, that is, £<°' d 0) 



= E. 



(0,0) 



is set to 1. Panels (a) and (b) are density plots of E^ d ' d and E^ , respectively, in the 
(k,l) space; darker regions correspond to lower entanglement. The thin lines are the isolines (lines of equal entanglement), 
while the thick straight lines correspond to k — I (dashed), k + I = 10 (dashed-dotted) , and I — 4 (dotted). The entanglement 
of the latter three lines is plotted in (c), (d), and (e), respectively, with dark- blue and light-red denoting the photon-added 
and photon-subtracted states, respectively. Even though k and I are physically discrete variables, we have taken them to be 
continuous variables by using the Gamma function as an analytic extension of the factorial function (this holds for the rest of 
the figures in the article). 



In order to prove that is a monotonically increasing allows us to write 

function of k, we proceed as follows. The Pascal identity „ „ n.,^ „ ,,n 

P^=\ 2 p^+{l-\ 2 W, (16) 

(n + k + 1\ _ (n + k\ (n + k\ 

\ k + 1 r HI k ' y ' where we set pit^ = for n < by definiteness. Now, 
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using the strict concavity of the function h[x] = — xlogx, 
we have 



n=0 



oo oo 

h[^]>X*^h[p<R)]+0--#)E h [ 



n=0 



n=0 



,(*0 



(17) 



for < A < 1. Since p^-i is equivalent to pn +1 ^ up to a 
shift to the right in the Fock basis, which does not change 

the entropy, i.e., 



n 

the expression (|17[) is simply equivalent to 



n=Q ' 



_(*+!) 
Pn-1 



(fc+1) 



(fc+1) 



(18) 



for < A < 1. Thus, we conclude that the entanglement 
can only increase with the number of photon additions 
or subtractions when acting on one mode only (before or 
after the two- mode squeezer). 



IV. ENTANGLEMENT PROPERTIES FOR 
TWO-MODE OPERATIONS 

The non-trivial expressions of p£ fe '^ and q^' 1 ^ have pre- 
vented us from doing an exhaustive analytical study of 
the entanglement properties of states @ and (0 when 
operating on both modes, that is, when k ^ and I =/= 0. 
Indeed the only interesting property that we have been 

able to prove analytically is that -E^dd = -^iub » * na ^ 1S > 
for symmetric operation (same number of operations on 
both modes k = I), additions and subtractions lead to 
the exact same entanglement, a result already noted in 
[39l ] in the k = I = 1 case. In order to prove this, note 
that by renaming the summation index as n = m + k, 
the subtracted state can be written as 



E 



An ( fc ) 



' ^i(fc + i,fc + i;i;A 2 ) 



|m, m) 



E \/Pm' k) \m,m), 



(19) 



which implies that iV'sub ) an< ^ IV'idd ) have the same 
Schmidt coefficients, hence the same entanglement. 

In the reminder of the section, we analyze the entangle- 
ment of the states ([2]) and (J5)) by numerically performing 
the sums (|7|) and (J8J), truncated at an upper limit that 
ensures that the distributions p^h' 1 ^ and g! are normal- 
ized up to an accuracy of 10 -10 . Our numerical results 
are plotted in Fig. |3] for A = 0.4 and Fig.Hfor A = 0.22. 

The main tendency we can deduce from our numerical 
analysis is that it is always better to perform addition 
rather than subtraction in order to increase the entan- 
glement, i.e., > ^W> • This * s c l ear ly visible when 
comparing the density plots in Figs. [3jt and [3Jd. Such a 
result could be linked to the fact that photon addition 
seems to increase the state's non-classicality faster than 
photon subtraction does [49|, [Bfj. Note, however, that for 



large squeezing the difference becomes less pronounced, 
i-e-^iub^^ddforA^l. 

It is also important to remark that the probabilities of 
success of the addition and subtraction schemes are dif- 
ferent [H, |3(| S3 , and hence, even though addition per- 
forms better for the same number of operations, it might 
be preferable to perform more subtractions to achieve a 
given entanglement, depending on the particular experi- 
mental scenario. 

For symmetric operation k = I, where addition and 
subtraction perform equally, the entanglement increases 
with the number of operations, i.e., E^ k+1 ' k+1 ^ > E^ k ' k \ 
This is explicitly shown in Fig. [3J:, where we plot the 
entanglement E^ k ' k ^ as a function of k. This behavior is 
in agreement with the studies performed in [4(J EH ■ 

We also observed that for a fixed number of opera- 
tions k + I, the entanglement increases when approach- 
ing the symmetric situation k = I. This is shown in 
Fig. [5Ji, where we plot E^ k,w ~ k ^ as a function of k both 
for addition (dark-blue line) and subtraction (light-red 
line). Note, however, that the shapes of the curves are 
rather different for addition and subtraction. Neverthe- 
less, we remark that the optimal enhancement is obtained 
when the same number of operations are applied on both 
modes, where both addition and subtraction give the 
same entanglement enhancement. 

By keeping the number of operations fixed on one 
mode, the entanglement is an increasing function of 
the number of additions on the other mode, that is, 



E (k+U) 

^add 



^add ■ Thus, by fixing I to some value, say 
^o, the entanglement increases as Alice adds more pho- 
tons; for Zo = this is exactly the analytical result that 
we found in Sec. III. The case Iq = 4 is illustrated in 
Fig. [3^. The situation is a bit different for photon sub- 
traction. While above some critical squeezing parameter 
A (depending on Iq), the entanglement is a monotoni- 
cally increasing function of the number of subtractions k 
for a fixed I = Iq, just as for additions, below this crit- 
ical squeezing it is not. This is made clear in Fig. 31 

(k I) 

where we plot ^ for a smaller squeezing parameter 
A = 0.22. Note Fig. 2J;, in particular, where we see that 
the entanglement decreases in some interval of k above 
the symmetric point k = I = 4 before going back to its 
normal increase. Otherwise, the behavior of entangle- 
ment at A = 0.22 as shown in Figs. 0k and|3}} is qualita- 
tively similar to what we observed at A = 0.4 in Fig. [3] 
The case of photon addition is also plotted in Figs. 2Jd 
and [4}: for comparison. 

Comparing Figs. [3] and 21 we also observe that the en- 
tanglement enhancement effected by photon addition and 
subtraction is greater, in relative terms, when the squeez- 
ing parameter A is low (remember that -E^dd'' and E^£ 
are normalized to the TMVS state in the figures). For 
example, at the symmetric point k = I — 4, the entangle- 
ment is enhanced by a factor about 3.7 with respect to 
the TMVS state at A = 0.4, while it reaches about 6.7 at 
A = 0.22. This can be interpreted as follows. For A —> 1, 
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FIG. 4: (Color online) Same as Fig. [3] but for A = 0.22. For conciseness, we only show as a density plot (a) as well 

as the entanglement along the lines k + I = 10 (b, dashed-dotted line) and I = 4 (c, dotted line), both for the photon-added 
(dark-blue) and photon-subtracted (light-red) states. Note that while the photon-added states have the same behavior as for 
A = 0.4, this is not the case for the photon-subtracted states, see details in the text. 



the entanglement of the TMSV state is already very large 
and its Schmidt coefficients approach a uniform, infinitely 
wide distribution; hence, photon addition or subtraction 
cannot improve much on this. 

Finally, it is worth comparing the photon- added and 

-subtracted states, IV'idd ) an< ^ l^sub )' f rom tnc point of 
view of the energy cost for generating the same amount 
of entanglement. For this, we define the entanglement 
energy-efficiency of the photon-added state as 



'/add 



E, 



(k,i) 

add 



where N. 



(k,l) 

add 



V^add 



a) a 



(^dd°/2) 



(20) 



add 



is the total mean 



photon number of the state, and the function g(x) — 
(x + 1) log(x + 1) — x\og(x) is the entanglement entropy 
of a TMSV state with 2x mean photon number (equal to 
the entropy of a thermal state of x photons). Of course, 
we use a similar definition for photon-subtracted states. 

Taking into account that the TMSV state provides the 
highest entanglement for a given average photon num- 
ber, the entanglement energy-efficiency as defined here is 
equal to one for a TMSV state and less than one other- 
wise. In other words, the efficiency quantifies the degree 
to which the state's energy is optimally deployed in cre- 
ating entanglement. 

In Fig. [5J we plot the entanglement energy-efficiency 
as a function of the number of operations (fc, Z), for two 



values of A. Note that, even though photon addition leads 
to a larger entanglement amplification in absolute terms 
as shown before, the results shown in Fig. [3] tell us that 
photon subtraction is more efficient (in general) from the 
energy cost point of view. 



V. NON-GAUSSIANITY OF THE 
PHOTON-ADDED AND -SUBTRACTED STATES 

In this section, we evaluate the non-Gaussianity of the 
photon-added and -subtracted states that we have intro- 
duced in the previous sections, and investigate the possi- 
ble link with their entanglement properties. In a nutshell, 
we reach the conclusion that photon addition leads to a 
faster degaussification of the TMSV state than photon 
subtraction, which is reminiscent to the behavior of en- 
tanglement, but nevertheless the level of entanglement 
found in these states seems to have no direct relation to 
their non-Gaussianity. 

We use here the non-Gaussianity measure of a state 
that was introduced in (5l| . This measure has already 
been used in a similar context, for example in [52j ■ There, 
after a numerical analysis based on this measure, it was 
conjectured that, at least for the class of photon-number 
entangled states (to which the states included in this 
work belong), the entanglement of Gaussian states is 
more robust against lossy channel with thermal added 
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FIG. 5: (Color online) Entanglement energy-efficiency rf k ' l > of the photon-added (dark-blue curves) and photon-subtracted 
(light-red curves) states as a function of the number of operations k and I for A = 0.4 (upper row) and A = 0.22 (lower row). 
As we did for the entanglement, we show the evolution of the efficiency along the k = I (a,d), k + I — 10 (b,e), and I — 4 (c,f) 
lines. 



noise than that of non-Gaussian states. Note, however, 
that this conjecture was recently proved wrong by show- 
ing that it does not hold when a different entanglement 
measure (negativity under partial transposition) is used 
[53[ | , and that the entanglement of the NOON states and of 
a simple class of photon-number entangled states survive 
longer in a thermal environment than the entanglement 
of any Gaussian state [54[ . 

Let us first explain how this non-Gaussianity measure 
G[p] works for a general state p 51] . The idea is to eval- 
uate the statistical distinguishability between p and the 
Gaussian state pq having the same first and second mo- 
ments, which can be done by using the quantum relative 
entropy. Thus, we define the non-Gaussianity of state p 



states have a covariance matrix of the type 
c (k,i) = 



r a (M) t 


7 (M) <J Z 1 


_ 7 (^) a z 





(22) 



where we have defined the 2x2 matrices 1 = diag(l, 1) 
and a z = diag(l, — 1). In the case of the photon-added 

states IV'adcP), the parameters are given by 



= l + 2fc + 2^>^>, 

00 

^ = 1 + 2Z + 2 £ »#.«>, 



(23) 



G[p] = S[p\\p G ] = tr{p(logp-log/3 G )} 
= S[p G ] -S[p\. 



(21) 



where the last equality follows from the fact that p and 
Pq have the same first and second moments. Here, the 

states l^idd ) an d IV'sub) wc wor k with are pure, and 
hence their non-Gaussianity is simply the entropy of the 
corresponding Gaussian state, i.e., G[p] = S[p G }- 

Now, let us define the vector operator f = 
(x a ,p a , Xb,Pb) built on the quadrature operators x a = 

+ a and p a = i(a' — a), and similarly for the mode 
B. It is fairly simple to check that the states | V'ad'd ) 
and iV'sub ) n ave all zero mean, that is, (f) = 0. The 
elements of the covariance matrix are then evaluated 
as Gji — (fjfi + f;fj)/2, and it is straightforward to 
check that both the photon-subtracted and photon-added 



/add 



= 2 J2 V(" + k + l)(n + I + l)^''^ 



n=0 



while, for the photon-subtracted states IV'sub )' one nas 



(k,i) 



-sub 



(k,l) 
7sub 



l-2fc + 2 

?i-max(fc,/) 



(24) 



1-2/ + 2 J2 n 1n' l \ 

n— max(fe,/) 



n— max(/c,/) 



While we have not been able to evaluate the sums in the 
7's analytically, the sums in the a's and (3's have the 
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FIG. 6: (Color online) Non-Gaussianity (5' fe ''' of the photon-added (dark-blue curves) and photon-subtracted (light-red curves) 
states as a function of the number of operations k and I for A = 0.4 (upper row) and A = 0.22 (lower row). As we did for 
the entanglement and the entanglement energy-efficiency, we show the evolution of the non-Gaussianity along the k — I (a,d), 
k + l — 10 (b,e), and I = 4 (c,f) lines. Note that, in contrast with the entanglement and the entanglement energy-efficiency, the 
non-Gaussianity has the same qualitative behavior for both A values. 



following closed expressions: 



n=0 



np 



(k.i 



) = (l + fc)(l + /)A 



;2 f 1 (fc + 2,/ + 2;2;A 2 ) 
2 F 1 (k + 1,1 + 1;1;A 2 )' 

(25) 



5>?^> = k + {i + k)\^ k 

n=k \l) 

2 F 1 (k + 2,k + 2;k-l + 2;X 2 ) 
* 2 Fi{k + l 1 k + l-k-l + l ] \*Y 

where in the second expression we have assumed k > I 
(once again, exactly the same expression but interchang- 
ing k and I holds for k < I). 

The two-mode covariance matrix (|22[) is in normal 
form Q, from which the entropy of the Gaussian state 
(and thus the non-Gaussianity of the photon-added or 
-subtracted state) can be directly evaluated as 



r (fc,Z)l 



where 



, , z + 1, z+1 z — 1 z — l 
#0) = — — log — — log — — . 



(26) 



(27) 



In Fig. |51 we plot the non-Gaussianity G^ k '^ for an ar- 
bitrary number of operations (fc, I) with A = 0.4 (upper 
row) and A = 0.22 (lower row). In analogy with the be- 
havior of entanglement, we observe that photon addition 
leads to a faster degradation of the Gaussianity of the 
TMSV state than photon subtraction. In other words, 
iV'id'd ) (dark-blue lines) is more non-Gaussian than the 
photon subtracted state l^suiP) (light-red lines). This is 
clear for example in Figs.[5fi,d where we plot the increase 
of & k ' k ^ with k for symmetric operations k = I. Extrap- 
olating from the behavior of entanglement, one would 
be tempted to predict that the non-Gaussianity £?'*>0 is 
maximum for symmetric operations k = I for both addi- 
tion and subtraction. Interestingly, the behavior of 
is radically different. While for a fixed number k + I of 
photon additions (dark-blue lines), its maximum is in- 
deed reached for k = I, for a fixed number k + l of photon 
subtractions (light-red lines), the non-Gaussianity is ac- 
tually minimum for k = I, see Figs. [5}D,e. We observe 
a similar anomaly in Figs. [BJ;,f, where we plot the non- 
Gaussianity as a function of k for a fixed I. Thus, the non- 
Gaussianity of photon-subtracted states exhibits, in some 
situations, a very different qualitative behavior from that 
of its entanglement. 



and where 



,(M) 




(k,l)2 _|_ 



(28) 



are the symplcctic eigenvalues Q of the covariance ma 
trix 



VI. CONCLUSIONS 

We have studied how local photon addition and sub- 
traction affect the entanglement and Gaussianity of the 
two-mode squeezed vacuum state. This subject has be- 
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come of recent interest, especially since these fundamen- 
tal heralded non-Gaussian operations have become avail- 
able in the laboratory. 

First, we have analytically shown that the entangle- 
ment grows with the number of photon additions or sub- 
tractions when only one of the parties performs the op- 
erations. Wc have then numerically analyzed the case in 
which both parties add or subtract photons; although ad- 
dition and subtraction lead to the same entanglement en- 
hancement when both parties perform the same number 
of operations, photon addition beats photon subtraction 
in general. 

We have also analyzed the efficiency with which the en- 
ergy in photon-added and photon-subtracted states gen- 
erates entanglement, showing that, in general, this can 
be close to perfect for photon subtraction, but not for 
photon addition. 

Finally, we have numerically studied the degaussi- 
fication of the two-mode squeezed vacuum state that 
is effected by photon addition or subtraction, show- 
ing that photon addition degrades the Gaussianity of 
the state faster than photon subtraction. Observing 
that the entanglement and no n- Gaussianity of photon- 
subtracted states have a radically different behavior, wc 
conclude that the relation between entanglement and 
non-Gaussianity is not as simple as previously assumed. 

Future directions might include analyzing how photon 
addition and subtraction affect the entanglement of more 



general (possibly mixed) states, such as the TMSV state 
degraded by losses in the channels through which the 
entangled modes are sent to Alice and Bob, which is a 
typical initial state of many distillation or concentration 
protocols {42j . 
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